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This  is  basically  a theoretical  paper  moti- 
vated by  a problem  from  seismic  data  processing 
in  oil  exploration.  We  develop  a Kalman  filtering 
approach  to  obtaining  optimal  smoothed  estimates 
of  the  so-called  reflection  coefficient  sequence. 
This  sequence  contains  important  information 
about  subsurface  geometry.  Our  theoretical  pro- 
blem is  one  of  estimating  white  plant  noise  for 
the  systems  described  in  Eqs.  (1)  and  (2).  By 
means  of  the  equations  which  are  derived  herein, 
it  is  possible  to  compute  fixed  - interval,  fixed- 
point,  or  fixed-lag  optimal  smoothed  estimates  of 
the  reflection  coefficient  sequence,  as  well  as 
respective  error  covariance -matrix  information. 
Our  optimal  estimators  are  compared  with  an  ad 
hoc  "prediction  error  filter,  " which  has  recently 
been  reported  on  in  the  geophysics  literature. 

We  show  that  one  of  our  estimators  performs  at 
least  as  well  as,  and,  in  most  cases  better  than, 
the  prediction  error  filter. 


s 
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determined  by  the  material  properties  of  the  me- 
dium. Whenever  such  a wave  encounters  a change 
in  acoustic  impedance  (acoustic  impedance  is  the 
product  of  rock  density  and  rock  propogation 
velocity),  a certain  fraction  of  the  incident  wave 
is  reflected  upwards.  Seismic  detectors  situated 
at  the  earth's  surface  record  the  continual  motion 
of  the  earth  under  the  impact  of  seismic  waves 
impinging  from  below.  This  recording  is  per- 
formed digitally  at  a fixed  sampling  increment. 

The  resultant  set  of  discrete  observations  is 
called  a . . . 'seismic  trace'  and  constitutes  a 
sample  of  a time  series." 

In  1954,  Enders  A.  Robinson  (4  and  5)  pro- 
posed a convolution  summation  model  to  describe 
the  signal  received  by  a seismic  sensor.  We 
write  this  model,  as 


z(k)  = VR(k)  + n(k) 

k 

v (k)  = r M(j)vV-j) 
R j=i 


(3] 


(4) 


INTRODUCTION 


This  is  basically  a theoretical  paper  moti- 
vated by  a problem  from  seismic  data  processing 
in  oil  exploration.  The  theoretical  problem  is 
one  of  estimating  white  plant  noise  for  models 
J.  and  : 

x(k  + 1)  - ixjk)  + jj(k) 
zjk)  = Hx(k)  + v(k) 

x(k  + 1)  = $x(k)  + Twlk) 
zjk)  - Hx(k)  f v(k) 

The  rest  of  this  section  briefly  describes  the 
seismic  data  processing  problem  and  its  relation 
to  this  theoretical  problem. 

One  of  the  most  useful  methods  to  predict 
the  presence  of  oil  is  reflection  seismology.  On 
land,  for  example,  an  explosive  is  detonated  be- 
low the  earth's  surface,  imparting  a pulse  of 
energy  into  the  earth.  * As  Wood  and  Treitel  (L 
state.  ''This  source  pulse  ...  is  split  into  a large 
number  of  waves  traveling  along  various  paths 


/jppy  »y*M  vmc  pi  p>r;  ;• 


where  V^(k)  is  the  noise  free  seismic  trace;  n(k) 
is  "measurement"  noise  which  accounts  for  physi- 
cal effects  not  explained  by  V^(k),  as  well  as 
sensor  noise;  k is  short  for  time  t,  ; V (i), 
i = 0, 1,  2,  . . . , I,  is  a sequence  associated  with  the 
basic  seismic  wavelet  (2^  and  5J;  and  u(j), 

j = 1,  2 is  the  reflection  coefficient  sequence. 

This  convolution  summation  model  can  be  derived 
from  physical  principles  and  some  simplifying 
assumptions,  such  as;  normal  incidence  (i.  e., 
horizontal  layering),  each  layer  is  homogeneous 
and  isotropic,  small  strains,  and,  pressure  and 
velocity  (or  displacement)  satisfy  a one-dimension- 
al wave  equation.  ^ Signal  V K ) , which  is  record- 
ed at  the  earth's  surface,  is  a superposition  of 


The  reader  interested  in  elements  of  the  seismic 
prospecting  method  and  the  seismic  reflection 
technique  should  see  Anstey  (2^),  Chs.  1 and  3, 
and  Dix  ( 3_). 

^ It  is  useful,  for  conceptual  purposes,  to  draw 
the  analogy  between  a layered  earth  system  and  a 
sequence  of  connected  lossless  transmission  lines. 
At  the  transmission  line  connections  (discontinui- 
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wavelets  which  are  reflet  ted  from  the  interfaces 
of  subsurface  layers.  The  u(j)  are  related  to 
interface  reflet  tion  and  transmission  coefficients. 

The  reflection  coefficient  sequence  often  has 
peaks  at  subsurface  interfaces,  and  is  often 
assumed  to  be  random  and  uncorrelated  ( 5 ).  It 
can  be  represented  graphically  as  a sequence  of 
''spikes”  which  occur  at  k = I,  2, ...  (Figure  1).  We 
would  like  to  learn  as  much  as  possible  about  the 
peak  values  of  the  p(j)  sequence,  including  ampli- 
tude and  time  of  occurrence,  so  that  information 
can  be  assembled  about  the  subsurface  geometry. 

If  noise  was  not  present,  and  we  had  numer- 
ical values  for  sequences  VR(i)  and  V + (i),  then 
the  reflection  coefficient  sequence  could  be  recon- 
structed by  passing  Vn(i)  through  an  "inverse 
filter.”  Let  VR(z),  V + (z),  and  p(z)  denote  the 
z-transforms  of  sequences  VR(i),  V + ( i) , and  p(i). 
Then,  Eq.  ( 2 ) can  also  be  written,  in  terms  of 
these  z-transforms,  as  VR(z)  =p(z)V+(z),  from 
which  we  see  that  n(z)  = VR(z)(l/ V + (z));  hence, 
when  VR(z)  is  operated  upon  by  the  "inverse 
filter"  I/V+(z),  we  obtain  p(z). 

In  the  noisy  situation,  it  is  common  practice 
to  extract  the  reflection  coefficient  sequence  using 
a technique  which  was  developed  by  Robinson  (4 
and  S_)  known  as  predictive  deconvolution.  The 
major  component  of  this  technique  is  a digital 
Wiener  filter,  which  requires  correlation  function 
information.  The  following  modeling  assumptions 
are  associated  with  Robinson's  work  (_1):  (1)  the 
layered  earth  is  a linear  system,  (2)  the  basic 
seismic  wavelet  is  minimum  delay  (i.e.,  minimum 
phase),  and  ( 3 ) the  reflection  coefficient  sequence 
is  random  and  uncorrelated. 

Because  the  Wiener  filtering  approach  to 
predictive  deconvolution  is  limited  by  these  model- 
ing assumptions,  which  may  not  always  be  valid, 
it  is  important  to  look  for  alternative  approaches 
which  will  permit  more  flexible  modeling  assump- 
tions. Kalman  filtering  is  such  an  alternative, 
and  has  recently  begun  to  receive  attention  from 
geophysicists  (7-9J.  In  this  paper  we  shall 
develop  a Kalman  filtering  approach  to  obtaining 
optimal  estimates  of  the  reflection  coefficient 
sequence. 

In  order  to  use  a Kalman  filter,  a state 
space  representation  is  required,  in  contrast  to 
the  convolution  summation  representation  in  Eq. 

(4^),  which  is  the  basis  for  use  of  a Wiener  filter. 
From  linear  system  theory  (10,  for  example),  it 
is  well  known  that  the  output,  y(k),  of  a linear, 
discrete-time,  time-invariant,  causal  system, 
whose  input  r(k)  is  zero  prior  to  time  zero,  is 

k 

y(k)  = T.  r(j)  h(k-j)  ( 5 ) 

j=0 

where  h(i)  is  the  unit  function  response  of  the 
system.  Comparing  Eqs.(5)  and  (4),  we  are  led 

2(contd) 

ties),  waves  are  reflected  and  transmitted,  just 
as  they  are  at  the  interfaces  between  earth  layers. 


to  tiie  following  important  system  interpretation  for 
the  seismic  tra.  e model;  signal  Vq(k)  i an  be- 
thought of  as  the  output  of  a linear  time-invariant 
system  whose  unit  response  is  y*(i)  and  input 
sequence  is  the  reflection  coefficient  sequence, 

-;(i).  (Figure  2), 

We  shall  assume  that  a state  space  represen- 
tation of  Eqs.  ( 3 ) and  (4)  is  our  starting  point  for 
predictive  deconvolution.  How  to  obtain  such  a 
representation  from  the  available  information  is 
not  the  subject  of  the  present  paper;  but,  is  an 
important  topic  that  will  be  discussed  in  the  con- 
text of  this  application  in  a future  paper. 

For  a single-input  single-output  system,  a 
most  efficient  state  space  realization  is  the  phase 
variable  canonical  form  (11): 

x(k  + 1)  = $xjk)  + bp(k)  ( 6 ) 


z(k)  = h x(k)  + n(k) 

where  x - col  (x  , x , 
li  = col  ( 1,  0,  ....  0), 


xn).  b = col(b1,b2 bn). 


and  Xj(k)  = VR(k).  Accepting  Eqs.  (6)  and  ( 7 ) as 
our  starting  point  for  predictive  deconvolution, 
instead  of  Eqs.  (3)  and  (4),  we  shall  study  the 
estimation  of  the  "white”  reflection  coefficient 
sequence,  M(k),  by  means  of  Kalman  filter  theory. 

Bayliss  and  Brigham  (2)  and  Crump  (jj)  have 
attempted  to  apply  kalman  filtering  to  predictive 
deconvolution;  however,  they  assume  p(k)  is  a 
colored  noise  process  and  proceed  along  a route 
which  is  quite  different  from  the  one  we  shall 
take.  Ott  and  Meder  (2)  have  attempted  to 
estimate  p(k)  via  Kalman  filtering;  however,  as 
we  point  out  in  Section  TV,  their  estimator  is  an 
ad  hoc  one,  and  does  not  perform  as  well  as  ours. 

We  shall  imbed  the  features  of  Eqs.  (6)  and 
(7)  in  a more  general  setting;  namely,  we  shall 
study  the  problems  of  estimating  white  plant  noise 
for  the  models  and  which  are  given  in  Eqs. 

(1)  and  (2)  respectively.  Optimal  estimators  for 
uj(k)  and  w(k)  are  developed  in  Section  II.  The 
relationship  of  zero  estimates  for  w(k)  to  zero 
Markov  parameters  is  examined  in  Section  III. 

We  relate  our  results  to  the  Ott  and  Meder  results 
in  Section  IV. 

OPTIMAL  WHITE-PLANT-NOISE  ESTIMATORS 
Introduction 

Let  us  direct  cur  attention,  for  the  time 
being,  at  Jy  in  Eq.  (1),  in  which  x6Rn,  weRn, 
£€Rm.  vf Rm,  and  ■i_(k)  and  v(k)  are  uncorrelated 
gaussian  white  sequences  with  covariances  Q and 
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K,  respectively.  Our  objective  is  to  obtain  a 
minimum-variant  e estimator  for  the  white-plant 
noise  sequence  j_(k). 

Let 

Z(j)  I z(  1 ).  z<  2 ) z(  j ) ) . ( 8 ) 


in  Z(j). 
theory 
x ( k t 1 1 k ) - 
Z(k  + 1)}. 
computed 


and  x(k)j)  denote  the  minimum  variance 
of  x(k)  which  uses  all  the  measurements 
From  a fundamental  result  in  estimation 
(12,  for  example),  it  is  well  known  that 
E]x(k  *l)|Z(k)]  and  x(k  + l|k  + l)  = E)x(k+l) 
and  that  x(k  + l|k)  and  Jf(k  f 1 1 k + 1)  can  be 
from  the  Kalman  filter,  whose  equations 

x(k  *■  1 |k)  = $x(  k| k ) 

P(k  t-  1 [ k)  = i P(  k | k )}’  + Q 

K(k  + 1)  = P(k  ^ llk)H’[HP(k  + 1 1 k)H*  + R]' 

£(k  + l|k)  = z(k  + l)  - Hx(k  h 1 1 k) 

x(k  + 1 1 k + 1)  = x(k  t l|k)f  K(k  + l)z(k  + 1 } k) 

P(k  + l|k  + 1)  = [I  - K(k  + I)H]  P(k  + 1 1 k) 


Observe,  from  Eq.  (1),  that  z(k  f 1)  depends 
on  x(k),  t£(k  - 1),  , u>(  1 ),  whereas  z(k)  depends 

only  on  x(k  - 1),  x(k  - 2),  ...»  a>(  1 );  hence  Z(k  + 1) 
depends  on  j,(k),  but  Z(k)  depends  at  most  on 
x(k  - l).  The  single-stage  predicted  estimate, 
x(kfl|k),  in  Eq.  (9)  can  be  obtained  directly 
from  Eq.  (1)  by  operating  on  both  sides  of  the 
state  equation  with  the  conditional  expectation 
operator  E|  . |Z(k)}.  The  term  E i^(K)  1 Z(k)}  = 

Ej  n(k)1  = 0 by  virtue  of  our  preceding  remark  and 
the  fact  that  r>(k)  is  a zero  mean  process. 

In  order  to  obtain  the  first  meaningful  mini- 
mum variance  estimate  of  ;r(k),  operate  on  both 
sides  of  the  state  equation  in  Eq.  (1)  with  the 
operator  Ef  - | Z(k  +2/}  ; 

E | x(k  tl)|  Z(k  cl)}  = i E|  x(k)  | Z(k  + 1)} 

+E|a)(k)|  Z(k  + 1)}  (15) 

The  last  term  in  Eq.  ( 15  ) will  not  be  zero  because 
Z(k  + 1)  depends  on  r_(k);  hence,  we  obtain  the 
following  estimator  for  w(k) : 

5(k|k  + l)  = x(k  * l|k  f 1)  - ix(k|k  + 1)  (16) 

Observe  that  this  is  a single-stage  smoothed  es- 
timate of  x(k)  which  requires  not  only  the  optimal 
filtered  estimate  of  x(k  + l),  but  also  the  optimal 
single-stage  smoothed  estimate  of  x(k),  for  its 
implementation. 

By  way  of  these  calculations  we  have  estab- 
lished the  fact  that  minimum  variance  estimates  of 
X(k)  will  be  optimal  smoothed  estimates.  In  the 
remaining  paragraphs  of  this  section  we  shall  ex- 
plore T(k|kfl)  and  its  estimation  error  in  more 
detail,  obtain  recursive  equations  for  an  optimal 
4 -stag**  smoothed  estimator,  £(k|k  + £),  and  its 
error  covariance  matrix,  and,  shall  obtain  compa- 
rable results  for  estimates  of  w(k)  in  Eq.  (2). 


Single  -Stage  Opt i rnal  Smoothed  Estimates  of  'ij(k) 
Mere  we  develop  a more  useful  equation  for 
T.(k|k<  ll  then  Eq.  and  develop  an  equation 

tor  the  error  covariance  matrix  «Ji  (klktl),  where 

•1) 

’i,A(k|k  t l)  - E|  [T( k | k ♦ 1)  - E|T(k  |k  i l))][same]  } ( 17  ) 
Theorem  1.  For  system  Jy 

i(k|k  t 1)  =QP'L(k  + l|k)K(k  + l)z(k  f l|k)  ( 18  ) 


*,t(k  |k  + 1)  r Q . QH’[HP(k  + 1 |k)H’  + R]  HQ  ( 19  ) 


Proof:  (a)  Derivation  of  Eq.  ( 18  ).  Meditch  ( 12  ) 
derives  the  following  equation  for  x(k|k  + l): 

x(k(k  + 1)  = x(k|k)  + A(k)K(k  + l)z(k  + l|k)  (20) 

where 

A(k)  - P(k|k)  t'p"\k  + l|k)  (21) 

Substitute  Eqs.  (20)  and  (13)  into  Eq.  (16),  making 
use  of  Eq.  (9),  to  show  that 

£(k|k  + l)  = Cl-tA(k)]K(k  + l)z(k+l|k)  (22) 

From  Eqs.  (21)  and  (10),  it  follows  that 


I - $A(k)  =QP"  (k  t l(k)  (23) 

Substitute  Eq.  ( 23  ) into  Eq.  ( 22  ) in  order  to 
obtain  the  desired  expression  for  i(k|ktl)  in  Eq. 
(18). 

(b)  Derivation  of  Eq,  (19).  Signal  ?(k  + l|k), 
the  innovations  process,  is  known  to  be  zero  mean 
and  white  (^3),  and, 

Eiz(k  + l|k)?,(k  + l|k)]  = HP(k  + llk)H'  +R.  ( 24) 

Subsequently,  E^(k(k  + 1)}  = 0,  which  means  that 
E|w(k|k  + 1)}  = Et^(k)]-E|i(k|k  + l)}  = 0,  and  therefore, 
that 


yklk  + 1)  = EiS(k|k  + l)5'(klk  + l)} 

= E l [^(k)  - i(k|k  + l)][iu(k)  - iL(k  |k  + 1)]') 

= Q - E | w(k )fi' (k | k + 1 ) ) -Eii(k|k  + l)a<'(k)) 

+ E{T(k|k  + l)5'(k|k +1))  ( 2 

From  Eqs.  (18),  (11),  and  (24),  it  follows  that 

E|ii(k|k  + l)®’(k|k  + l)]  = QH'[HP(k  + l|k)H,tR]‘iHQ 


Next,  we  shall  show  that 


E|w(k)  z'(k  + l|k)}  Qll'  . 


J 
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Recall,  that 

£(k  + l|k)  H*(k  f l|k)  + v(k  + 1)  ( l\ 

where  ?(k  + 1 1 k)  x(k)  - x(k  < 1 1 k).  Mcditch  ( ) 

derives  the  following  equations  for  ?(ktl|k)  and 
x^(k|k)  [Eqs.  (5.63)  and  (5.64),  respectively  in 

Hi: 

x(k  + ]|k)  ^ i x(k|k)+(U(k>  (25 


£(k \K)  = [I  - K(k)H]4  x(k  - l|k  - 1) 

+ [i- K(k)H]w(k  - 1) -K(k)v(k)  (30) 

Substitute  Eq.  (29)  into  Eq.  (28),  to  show  that 

z(k  + l|k)  = H$x'(k|k)  + H£(k)  + v(k  + 1)  ( 31 ) 

Observe,  from  Eq.  (30),  that  x(k|k)  does  not 
depend  on  w(k);  it  depends  at  most  on  £(k  - 1);  hence, 

E|jj(k)  z'(k  1 1 1 k) } = El£(k)x'(k|k)}(HJ)' 

+ QH1  + E|£(k)v'(k  + l)l 
= QH1  (32) 

which  is  Eq.  ( 27  ). 

Finally,  to  obtain  ijk|k  + l)  in  Eq.  (19),  sub- 
stitute Eqs.  (26),  (18),  and  (32)  into  Eq.  (25).  O 
Comment  1.  It  is  interesting  to  observe,  from 
Eq.  ( 18  ),  that  the  estimator  of  white  noise  process 
£(k),  3>(k|k  + l),  is  itself  a white  noise  process. 

This  is  especially  useful  for  the  predictive  decon- 
volution problem  that  was  described  in  Section  I, 
where  the  reflection  coefficient  sequence,  which  is 
to  be  estimated,  is  white.  Additionally,  we 
observe  from  Eq.  (19)  that  * (k|k  + l)£Q,  which 
means  that  we  can  expect  to  'So  better  by  using 
£(k|k  + l)  then  by  using  a zero  estimator  for  w(k). 

2-Stage  Optimal  Smoothed  Estimates  of  w(k) 

Here  we  generalize  the  results  of  the  pre- 
ceding paragraph,  obtaining  i(k|k  + i)  and  t (k|k  + 2). 
These  estimator  equations  are  useful  in  that  they 
permit  us  to  compute  fixed-interval,  fixed-point, 
or  fixed-lag  smoothed  estimates  of  the  plant  noise. 
We  define  i (k|k+2)  as  follows: 

♦ (k  | k + 2 ) - E ( f£(k  | k + 2)  - E | £(k  | k + 2)  ] ] [ s ame  ] 1 ) 


£(k  |k  + l)  - &(k)  - ft(k|k  + 2)  ( 

Theorem  2.  For  system  Jy 

£(k|k  + 2)  - £(k|k  + 2 - 1)  4 N(k|k  + 2)z(k  + i|k  + 2 -1) 

a ( 

where  2 = 1,  2, . . . , 5>(k  | k)  S 0 , 

N(k|k  + l)=QP‘*(k  + l|k)K(k  + l)  ( 


k t Jt -1 

N(k |k  + 2)  - QP*  (k  + 1 1 k)  n A(i)K(k  + 2) 
i=k  + l 
.3 


A(i)=P(i|i)*'P'1(i  + l|t) 


Additionally, 

i^lklk  + = k)k  + 2 - 1)  - N(k  |k  + i)[HP(k  + 2jk  + 2 -1)H* 

+ R]N'(k  | k + 2)  (39) 

where  ^(k  | k ) ^ Q. 

Comment  2.  Assuming  the  truth  of  Eq.  (35), 
then  another  representation  for  £(k|k  + 2),  which  is 
nonrecursive,  is 

2 

ui(k  |k  + 2)  = £ N(k|k  + j)z(k  + j|k  + j -1)  (40) 

3=1 

This  result  is  obtained  by  iterating  Eq.  ( 35  ) on 

2 . o 

Comment  3.  Suppose  we  have  a complete  set  of 
measurements,  £(1),  £(2),...,  £(L),  available. 

Then  we  can  use  all  of  these  measurements  to 
obtain  fixed-interval  estimates  of  £(  0 ),  tu(  1 ), . . . , 
w(L-l).  In  Equation  (40),  set  k + 2 = L and  let 
2 = 1,  2, . . . , L to  obtain  the  following  fixed-interval 
estimator  oi  plant  noise: 

2 

Ii(L-2 1 L)  = F N(L-2|L-2+j)T(L-2+j|L-2  + j-l)  (41) 

j=l 

A fixed-point  estimator  of  £(k)  is  £(k|k  + 2) 
where  k is  fixed  and  2 is  varied.  Our  formulas 
in  Theorem  2 are  in  fixed-point  format,  and  can 
be  used  to  enhance  estimates  at  specific  values  of 
k. 

A fixed -lag  estimator  of  £(k)  is  i(k|k  + 2) 
where  2 is  fixed  and  k is  varied.  As  such, 

Jj(k|k  + 2)  utilizes  the  window  of  measurements, 

£(k  + 1),  £(k  + 2),  . . . , £(k  + i).  For  example,  £(l|l+2) 
utilizes  measurements  £(2),  £(3),...,  £(2  + 1), 
whereas  i(2|2  + 2)  utilizes  the  measurements  £(3), 
£(4),...,  £(2+2).  A fixed -lag  estimate  of  w(k) 
would  be  useful  in  those  situations  where  we  decide 
to  use  only  a fixed  length  of  data  to  obtain  the 
optimal  estimate,  o 

Proof  of  Theorem  2;  (a)  Derivation  of  Eq.  (35). 

Let 

Z(k  + 2 -1)=  |£(1),  £(2) £(k  + 2-1)]  (42) 


Z(k  + 2)=  |£(1),  £(2) £(k  + 2)l 


n A(i)  = A(1  ) A(  2). . . A(I-1  )A(  I). 
i=l 


Op(*rat«  on  both  sides  of  the  state  equation  in  Kq. 
(I)  with  operators  E | • | Z(k  ♦ Z - 1) } and  E|*|Z(kt£)} 
to  obtain  the  follow  ing  results; 


r(k|k  k i-1)  x(k  k 1 1 k k 4-1)  - 4x(k | k k i -1) 

and 

7i(k  | k k X ) - x(k  k 1 1 k k i ) - t!c(k  | k k Z ) 


(44  ) 


(45) 


Meditch  (12)  derives  the  following  equation  [Eq. 
(6.  60)  in  12]  for  x(k  | j),  j - k k 1,  kk  2,  : 


x(k  | j)  - x(k|  j - 1)  k M(k  |j)  z(j|  j - 1) 
where 


(46) 


M(k|j) 


= n A ( i J 


K(j) 


(47) 


We  u9e  Eq.  (46)  to  obtain  expressions  for 
x(k|k  + i)  and  j£(k  f 1 1 k + £),  as  follows; 


x(k  | k k X)  - x(k)j) 
and 

x(k  k 1 1 k k X)  = x(k | j) 


(48) 


j = kkX 


(49) 


k = k k 1,  j = k k X 


Substitute  Eqs.  (48)  and  (49),  using  Eqs.  (46), 
(47),  and  (44),  into  Eq.  (45)  to  show  that 


kkX-1 

7i(k  |k  k X)  -l(k|k  k i-1)  k [i  - $A(k)]  U A(i) 

Uk  kl 


*K(k  kX)?(kk  X|k  k X - 1) 


(51) 


♦^(k  | k k X)  = K|5[(klk  -k  X)Tj'(k|k  k X)} 
where,  from  Eqs.  ( 34  ) and  ( 35  ), 
jT(k  1 k k X ) = T(k  | k k X - 1 ) - N(k|k  k X)£(k  k / | k + X-  1)  (52) 


Substitute  Eq.  (52)  into  Eq.  (51),  using  Eq.  (24) 
for  k = k f f,  to  obtain  the  following  expression  for 
* (k  | k k X): 


k R]N(k  |k  k X)  (53) 

-N(k|k  k X)E|z(k  k X|k  k X-1)TT' (k  |k  k X-  1)} 

- E)u<(k|k  k X - l)z'(k  k X | k k X-  1) )N*(k  | k t X) 


Hec  ause  z is  a white  sequence. 


Ejz(k  I j |k  kj-l)7'(k  I X|k  kX-l)J  - 0 


for  all  j tl  \ hence,  using  Bjqs.  (40)  and  (34),  it  is 
straightforward  to  show  that 


EjT(k|k  k X-  l)z'(k  k X|k  t X -1 ) } 

= E< Jj(k)  z'(k  k x|k  k X-  1)1 


(54  ) 


The  quantity  on  the  right-hand  side  of  Eq.  (54)  is 
evaluated  in  Appendix  A,  as  (X  | X-  1) : 


♦ ^ </  | X - 1)  = N(k  |k  k X)rHP(k  k X |k  k X - l)H'k  R]  ( 55  ) 


Apply  Eqs.  (54)  and  (55)  to  Eq.  (53)  to  obtain  our 
results  in  Eq.  (39). 


Optimal  Smoothed  Estimates  of  w(k) 

We  now  direct  our  attention  at  estimation  of 
w(k)  in  ^2,  Eq.  (2),  when  wf  Rq,  qSn.  Signal 
w(k)  is  white  with  covariance  Qj  . 

Let  w(k|kkX)  denote  the  minimum- variance 
X-stage  smoothed  estimator  of  w(k),  and  ^(klkkX) 
its  error  covariance  matrix. 


*Jk|kkX)  = E|[w(k|kkX)-Ejw(k|kkX)}][  same] ' ] 


(56) 


Theorem  3.  For  System  , 

w(k)k  k X)  = w(k|k  k X - 1)  k N^(k  |k  k X)  |f(k  k X |k  k X - 1) 

(57) 


where  X=l,  2 w(k|k)  = £. 


( 50  ) Nw(k|kkl)  = Qir,P'1(kkl|k)K(k  + l) 


( 58  ) 


Apply  Eq.  (23)  to  Eq.  (50),  and  define  matrix 
N(k|kkX)  as  in  Eq.  (37),  to  obtain  our  Eq.  (35). 
We  define  N(k|k)  as  in  Eq.  (36)  so  as  to  be  able 
to  obtain  5(k|kkl)  from  Eq.  (35)  that  is  consistent 
with  our  earlier  results  in  Theorem  1. 

(b)  Derivation  of  Eq.  ( 39  ).  From  Eq.  (40), 
we  see  that  E f i(k  [ k k X)  ] = 0 ; subsequently, 

E |T(k  | k k X)}  = 0 , and 


N (k|kkX)=Qr'p1(kkl|k)  n A(i)  K(k  k X)  (59) 
w Uk+ 1 


kkX-1 


and  A(i)  is  defined  in  Eq.  ( 38 ).  Additionally, 
yk|kkX)  = 4w(k|kkX-l)  -Nw(k|k  + X)[HP(kkX|kkX-l)H' 
kR]N^(k|kkX)  (60) 

where  *w<k|k)*  Qj . 


Proof : System  J , in  Eq.  ( 1 ),  can  be  made 
equivalent  to  system  , in  Eq.  ( 2 ),  by  setting 
is(k)  = rw(k),  in  which  case 


+ (k  | k k X)  - * (k|k  k X- 1) + N(k|k  k XJTHPIk  k X|kk  X-1»h' 

•JJ  'J> 


w(k  |k  k X)  =Tw(k|k  k X) 
and 


(61  ) 


q = rQ,r 


(62) 


All  of  the  results  in  Theorem  3 follow  upon  sub- 
stitution of  Eqs.  (61)  and  (62)  into  the  respective 
Theorem  2 equations.  O 


646 


i 


i 


.A.. 


!• 


C onm.ci.t In  the  reflet  tion  Seismology  problem 

described  in  Section  1 Tsee  Fqs.  ( 6 ) and  ( 7 )J, 

w(k)  is  a Si  alar the  reflection  coeffi*  it  nt 

sequence,  hence,  the  results  in  Theorem  3 

are  the  ones  which  are  applicable  to  that  problem. 

ZLRO  MARKOV  FA  RA  M ETKRS  CASK 


( urollary  1.  Under  the  conditions  of  Theorem  4, 

V ( k 1 1<  t / ) Q,  for  i 1,<? i ( »*H  i 

w J 

The  proof  of  this  simple  result  follows 
directly  from  the  fat  t that,  bec  ause  of  Eq.  (66), 
w(k  |k  + l)  = w(k)  for  l 1,2,...,  j. 


Ott  and  Meder  (^)  consider  an  example  which 
they  state  ...  is  of  special  interest  for  seismic 
exploration  . . . ".  It  is  a one  dimensional  damped 
harmonic  oscillator  which  is  excited  by  impulses 
of  random  intensity  in  random  time  instances. 

For  their  example  , 


Corollary  2.  Under  the  conditions  of  Theorem  4, 

*w(k|k+j+l)  Q1-Q)(lKJr)’[HP(k  -t  j 4 1 J k + j)H ' 

tR]'1(HtV)Q1  (69) 

Proof:  From  Eqs.  (CO)  and  (C8),  it  follows  that 

*w(k|k  + j + 1)  = Q - Nw(k|k  + j + l)[HP(k  + j + l|k  + j)H' 

+ R]N^(k|k  + j + l)  (70) 

An  equation  for  N (k|k+j  + l)  is  developed  in 
Appendix  B,  Eq.  ^B-10).  Substitute  that  equation 
into  Eq.  (70)  to  obtain  the  desired  result  in  Eq. 
(69). 


If  one  substitutes  the  i,  F,  and  H matrices  for  this 
example  into  Eqs.  (57)  and  (58),  for  4=1,  he 
finds  that  £(k | k + 1)  = 0.  Additionally,  u(k|k  + 4)*0 
for  t » 2,  This  result  suggests  that  perhaps  the 
structures  of  },  r,  and  H establish  a first  value 
of  4 for  which  w(k|k+4)*0. 

Let  us  take  a closer  look  at  N (k|k  + l)  in 
Eq.  (56).  Substitute  Eq.  (11)  into  £q.  (58),  to 
shew  that 

Nw(k|k-  1)  =Q1(HF)'[HP(k  tl|k)H'  iR]’1  (65) 

Observe  that  N ( k | k + 1 ) depends  on  the  Markov 
parameter  HF  5 ).  In  the  Ott/ Meder  example 
HF  = 0 . 

We  generalize  the  preceding  observation  in 
the  following  theorem,  the  proof  of  which  is  given 
in  Appendix  B. 

Theorem  4.  For  system  J if  HT  = HiF  = ...  = 
H^‘V-0,  and  mV  < 0,  then 

w(k|kti)--0  for  1=1,2 j (66) 

and 

w(k  | k + j + 1)  * 0 . (67) 

This  theorem  states  that  if  the  first  j Markov 
parameters  are  zero,  then  the  first  non-zero 
estimator  of  w(k ) is  the  one  which  looks  j + 1 points 
into  the  future. 

Some  consequences  of  Theorem  4 are  given 
in  the  following  corollaries  to  that  theorem. 


Equation  (63)  is  obtained  by  discretizing  the 
equation  x(t )*er x(t)  + 8x(t)  = <u  u(t),  assuming  that 
t,  ,-t,  T 'S  small.  *" 

k = l k 


AN  AD  HOC  ESTIMATOR 

Ott  and  Meder  (^)  define  a 'prediction  error 
filter1'  for  j,(k),  in  system  J , as  the  difference 
between  x(k+l|ktl)  and  x(k  + l|k).  In  this  section, 
we  shall  examine  their  estimator  and  compare  its 
performance  with  our  estimator  l(k | k + 1 ((Section  II). 
We  shall  also  study  the  consequences  of  a pre- 
diction error  filter  for  estimating  w(k)  in  system 
As  such,  this  section  serves  to  illustrate 
some  of  our  theoretical  results. 

Ad  Hoc  Estimator  of  x( k) 

Combine  Eqs.  (13)  and  (9),  to  show  that 

x(k  + 1 |k  + 1)  = tx(k  |k)  + K(k  + 1)  ?(k  + 1 |k)  ( 71 ) 

in  which  we  recall  that  z“(k  + l|k)  is  the  white 
innovations  process.  Comparing  the  structures  of 
Eq.  ( 71  ) and  the  state  equation  in  Eq.  ( 1 ),  we 
define  a "prediction  error  filter"  for  £(k),  denoted 
here  as  5^(k),  as 

ij(k)  = K(k  + 1)  ?(k  + 1 1 k)  (72) 

Observe  that  ii  (k)  is  a white  estimator  of  ju(k), 
and  that  l^(k)  = x(k  + 1 |k  t 1)  - x(k  + 1 1 k).  Observe, 
also,  that  i(k)  is  not  an  optimal  estimator  of  ^(K ) ; 
it  has  merely  been  defined  as  in  Eq.  (72).  As 
such,  we  refer  to  the  prediction  error  filter  as  an 
ad  hoc  estimator. 

We  shall  compare  I^(k)  with  i(k|k  + l).  Let 
■5  (k)  denote  the  error  covariance  matrix  for  the 

a^*  hoc  estimator  uij(k);  i.  e.  , 

♦^(k)*  E|  [^(k)  - E|^1(k)1][sameV]  (73) 

where 
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i (k) 


-1 


(kl 


E | ~ (k > 1 0 , so  that 

— 1 — 

*ll*k|  t:l  i|(k) ) 


— l' 


(7S  ) 


Substitute  Fqs.  (74)  and  (72)  into  Eq.  (75),  using 
Eqs.  ( 24),  ( 27  ),  and  (11)  to  show  that 


•J  (k)  • Q + PH'(HPH’+  R)'lHP  -Qh'(HPH'+  R)'‘hP 
1 

- PH  (HPH  + R)  HQ  (76) 

in  which  P is  short  for  P(k  + l|k). 


Theorem  5.  For  system 


4x(k|k  +1)  * *ijul(k) 


(77) 


Tw^k)  = K(k  + l)z(k  + l|k) 


(79) 


w (k)  = (r'D’V'Ktk  + l)£(k  + l|k) 


(80; 


Let  t j(k)  denote  the  error  covariance  matrix 
for  wik),  defined  analogously  to  t'^fk)  in  Eq.  (73). 

Proceeding  as  we  did  for  the  calculation  of 
'll  (k)  in  Eq.  (76),  it  is  straightforward  to  show 
tftat 

♦wl(k)=Q1+(r’r)‘T,PH'(HPH'  + R)"1HPT(r,n‘l 


■ Q r'H'(HPH'+  R)'lHPV{T'?)'1 
- (T'r)'1r'PHl(HPH'+  R)"^H"Q. 


(81) 


*w(k|k^)Mwi(k) 


(82) 


Proof:  From  Eqs.  (76)  and  (19),  it  follows  that 
* ,(k)-t  (k|  k + 1)  = (P  - Q)  [h'(HPH'+  R)  "*H](P  - Q)  (78) 

JL)1  1! 

Matrix  R >0,  and  (HPH'+R)>0  for  the  inverse  of 
(HPH'+R),  which  is  used  to  calculate  K(k  tl),  to 
exist.  It  is  straightforward  to  show,  therefore, 
that  the  right-hand  side  of  Eq.  (78)  is  positive 
semi-difinite;  hence,  the  truth  of  Eq.  (77).  O 

The  significance  of  Eq.  ( 77  ) is  that  it  implies 
that  one  will  usually  obtain  better  performance  with 
our  single-stage  estimator,  i>(k|k  + l),  then  with  Ott 
and  Meder's  ad  hoc  predictic  >tror  filter;  and, 
that  one  cannot  do  worse  u .6  jur  estimator 
rather  than  the  prediction  error  filter. 

Ad  Hoc  Estimator  of  w(k) 

Now  compare  the  structures  of  Eq.  (71)  and 
the  state  equation  in  Eq.  (2),  to  obtain  the  follow- 
ing "correspondence1'  : 


In  this  equation,  w (k)  denotes  the  prediction  error 
filter  for  w(k).  Unfortunately,  Eq.  (79)  is  an 
overdetermined  system  of  equations;  i.  e.  , 
dim  Tw ,(k)  = n and  n>dim  wq(k) = q.  We  use  the 
pseudo-inverse  of  T,  (I'T)**’!'1,  to  obtain  the  follow- 
ing ad  hoc  prediction  error  filter  for  w(k) : 


The  proof  of  this  result  parallels  the  proof 
of  Theorem  5,  and  is  left  to  the  reader.  Once 
again,  we  see  that  we  can  expect  to  do  at  least  as 
good  using  our  single-stage  smoother,  w(k|k  + l), 
instead  of  the  prediction  error  filter,  Wj(k). 

Finally,  let  us  reconsider  the  harmonic 
oscillator  example  considered  by  Ott  and  Meder, 
whose  equations  are  given  in  Eqs.  (63)  and  (64). 
Recall,  from  our  discussions  in  Section  III,  that 
HP  = 0 for  their  system.  Under  this  condition, 

♦ (k)  in  Eq.  (81)  reduces  to 

t ,(k) ^q,+ (r,r)'1r,PH,(HPH,+ Ri^HprtrT)"1  (83) 

wl  1 

which  means  that,  for  their  example, 

*wl<k>aQl  (84) 

From  an  estimation  error-covariance  point  of 
view,  Eq.  ( 85  ) implies  that  the  prediction  error 
filter  performs  no  better,  and  usually  worse,  than 
the  zero  estimator  of  w(k)  I Why  bother  with  w (k), 
when  0 appears  to  be  a better  estimator  of  w(k)  y 
The  answer  to  this  question  is  related  to  the  seismic 
data  processing  problem  described  in  Section  I. 

We  are  interested  in  a time-series  for  the  reflec- 
tion coefficient  sequence.  By  examining  peak  values 
and  their  times  of  occurrence,  we  can  infer  some 
things  about  subsurface  geometry.  Even  a time- 
series  estimate  such  as  w (k)  will  have  peak  values 
(for  low  signal  to  noise  ratios)  that  make  w^(k) 
useful  [see  9 for  some  numerical  examples  that 
support  this  statement];  whereas,  a zero  estimate 
of  w(k)  does  not  contain  any  of  this  useful  informa- 
tion. 

CONCLUSIONS 

We  have  demonstrated  how  an  important 
application  can  be  interpreted,  from  a state  equa- 
tion point  of  view,  as  one  of  estimating  plant  noise. 
Our  emphasis  In  this  paper  has  been  on  the 
theoretical  development  of  white,  plant  noise  esti- 
mators. We  have  shown  that  such  estimators 
provide  us  with  smoothed  estimates  of  the  plant 
noise.  By  means  of  the  equations  which  are 
derived  in  the  main  body  of  this  paper,  it  is 
possible  to  obtain  fixed-interval,  fixed-point,  or 
fixed-lag  optimal  smoothed  estimates  of  white  plant 
noise,  as  well  as  respective  error  covariance- 
matrix  information. 

We  have  also  compared  our  optimal  estimators 
with  an  ad  hoc  "prediction  error  filter,"  and  have 
shown  that  better  performance  can  be  attained  with 
our  estimators. 

Much  work  remains  to  be  done,  including 
application  of  our  smoothers  to  real  seismic  data; 
study  of  computational  requirements  of  our 
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smoothers;  development  of  fast  realizations  of  our 
smoothing  equations;  and  the  extension  to  multi - 
input  multi-output  seismic  data  processing  problems. 
Recall,  also,  that  a major  feature  that  the  Kalman 
filter  has  to  offer  over  the  Wiener  filte  r is  it  can 
be  used  to  estimate  states  in  time -vary ing  systems. 
Additionally,  it  can  be  "extended"  to  estimate- 
states  in  nonlinear  systems.  The  extensions  of  the 
results  in  this  paper  to  estimation  of  the  reflection 
coefficient  sequence  for  time -varying  and  nonlinear 
source -earth  models  are  other  areas  for  additional 
work. 
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APPENDIX  A 

Calculation  of  E j2[(k  + 2 1 k + 2 - 1)  V (k)  ] 

For  notational  convenience,  let 


E|z(k+  i|k  + 2-  l)Jj'(k)  1 = yt|l-l)  ( A-l  ) 

Theorem.  For  system  J.  , 

(llO)  = HQ  (A-2) 

Z'U 

t (f  li-1)  = H4$*(k  + 2-l,k)»"1Q  (A-3) 

Tfuu 

where 

1 

4*(k  + 2 - 1,  k)  = II  |[I-K(k4j)H]*J  (A-4) 

j=2-l 

and  2 = 2,  3 Additionally, 


*'  (2|2-l)=N(k|k  + 2)[HP(k  + 2|k  + 2-l)H’+R]  ( A - 5 ) 

ZJJ 

where  N(k|k  + i)  is  defined  in  Eqs.  (36)  and  (37). 

Proof:  (a)  Derivation  of  Eq.  (A - 2) ( A =1  case). 

*»,  (1(0)  = E(?(k  + l|k)tj,'(k)]  ( A-6  ) 

The  truth  of  Eq.  (A-2)  now  follows  from  our 
developments  in  Eqs.  (27)  through  (32),  since 
il  (1 1 0)  is  simply  the  transpose  of  the  results  in 
Eq.  ( 27  )• 

(b)  Derivation  of  Eqs.  (A-3)  and  (A-4). 

Set  k = k + 2-  1 in  Eq.  (31)  to  show  that 

3f(k  + 2 |k  + 2-  1)  = Htx(k  + 2-l|k+i-l)  + H£(k  + 2-l) 

+ v(k  + 2);  ( A - 7 ) 

hence, 

(2|2-l)  = ME|x(k  + 2-l|k  + 2-l)w‘(k)}  (A-8) 

Z'U  — — 

since  E|w(k  + 2-  1) jj’ (k) J = 0 and  E |v(k  + 2)uj'(k)}  - 0. 
Next,  write  Eq.  ( 30  ) as 

x(k|k)  = **(k,k  - l)x(k  - 1 |k  - 1)  + t*(k,  k - Dt'^k  - 1) 

- K(k)  v(k)  ( A-9  ) 
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4 (k.  K - 1)  Cl-  K ( k ) 1 1 ] 4 (A  -10 

Thi'  solution  of  Eq.  (A-9)  (or  ‘x(k|k)  is  (Chapter 
1,  H) 

k 1 

5?(k|k)  - 4*(k,0)jf(0|0)  +£  4 *(k,  i)[  4 v(i,  i - 1)  4 uu(i-l) 
i = l 

- K(i)  v(i)3  (A  -11] 


4 *(k,  i)  - 4 *(k,  k-1)  4<'(k-i,  k-2) . . . 4*(iH,i)  (A -12> 

Set  k = k + 4-l  in  Eq.  (A -11)  and  observe  that  the 
only  term  which  depends  on  s(k)  is  the  one  for 
which  i = k > 1;  that  term  is 

4*(k  + 4- l.k  + 1)  4*(k  + 1.  k)  * _1^(k)  , 

which  can  also  be  written  as 

4 <‘(k  + 4 - 1,  k)  4 *w(k).  Based  on  these  results, 
evaluate  the  right-hand  side  of  Eq.  ( A - 8 ) to  show 
that 


(4|4-  1)  = Ht4*(k  + 4-l,k)  4 Q 


(A -13) 


which  is  Eq.  ( A - 3 ) . Substitute  Eq.  (A  - 10 ) into 
Eq.  ( A - 1 2 ) to  obtain  Eq.  (A-4). 

(c)  Derivation  of  Eq.  (A-5). 

Set  k = k + j in  Eq.  (14),  to  show  that 

I - K(k  + j)H  = P(k  + j |k  + j)P  (k  + j |k  + j - 1).  (A -14) 

Substitute  Eq.  (A  - 14 ) into  Eq.  (A-4),  take  the 
transpose  of  the  resulting  expression,  and  use  Eq. 
( 38  ) to  show  that 

4^(4 1 4-  1)  = QP‘\k  + l|k)  A(k  + 1)  A(k  + 2)  . . . A(k  + 4-  2) 
*P(k  + i- l|k  + 1-  1)  4'H'  (A -15) 

Next,  using  Eq.  (11)  in  which  k = k+4,  write 
N(k| k + 4)  as 

N(k  |k  + 4)  = QP'\k  + l|  k)  A(k  + 1)  A(k  + 2)  . . . A(k  + 4-  2) 
*P(k  + 4-l|k  + 4-l)4'H’[HP(k  + 4|k  + 4-l)H' 
+ R]*1  (A  -16) 

Comparing  Eqs.  (A  - 15 ) and  (A  -16),  we  see  that 


ami  (59): 

4 

w(k|k  I 4)  F N (k|k  < i)  £(k  t i|k  ti-1) 
i- 1 w 


I.cmma.  For  system  J , if  Hr  = H4T'=...  H4  T 

= 0,  then  N (k  I k t q)  0. 
w 

Proof:  From  Eq.  (50),  with  I=q, 

1 k + q-l 

N (k|k  + q)=Q  f P~  (k  + l|k)  A(m)K(k  + q)  (B-2) 

W m=k  + l 

The  following  equation  (12),  which  is  a well-known 
alternative  to  Eq.  (14),  is  used  in  our  proof: 

P_1(k  + l|k)  = P"*(k  + l|k+  1)  - h'r_1H  (B-3) 

Substitute  Eq.  (B-3)  into  Eq.  (B-2),  and  assume 
HT  = 0 , to  show  that 

k + q-l 

N (k|k+q)  = Q.r’P  (k  + l|k  + l)  n A(m)K(k+q) 

W m=k+l 

(B-4) 


, 1 k + q-l 

Nw(k|k  + q)  =Qjr  P"  (k+l|k+l)A(k+l)  H A(m)K(k  + q) 

m=k+2 

(B-5) 

Substitute  Eq.  (38),  for  i = k+l,  into  Eq.  (B-5)  to 
show  that 

k + q-l 

Nw(k|k  + q)=Q  rVP’  (k  + 2|k  + l)  n A(m)K(k  + q) 

m=k+2 

(B-6) 

Substitute  Eq.  (B-3),  for  k = k + l,  into  Eq.  (B-6), 
assume  H4T  = 0,  and,  substitute  Eq.  ( 38 ),  for 
i = k + 2,  into  the  resulting  expression,  to  show  that 
» 2 1 k + q-l 

N (k|k  + q)=Q  T (*  )'P"  (k+3|k  + 2)  H A(m)K(k+q) 
W m=k+3 

(B-7) 

Continue  this  development  assuming  H$^r  = 0,  then 
= 0,  • . • , and  then  F = 0,  to  show  that 

Nw(k  | k + q)  = Q jf'  (4q'1)'p'1(k  + q | k + q -1)  K(k  + q) 


Finally,  substitute  Eq.  (11),  for  k=k  + q,  into  Eq. 
(B-8)  to  show  that 


N(k|k  + 4)  = *:  (4|4-l)[HP(k  + 4|k+4-l)H’+R]'1.  . , q-1  . . . -1 

*JI  N (k  k + q)  = Q r (4q  l)  H tHP(k  + q k + q-l)H  +R] 

/A  171  W 1 


(A -17 ) 

from  which  Eq.  (A-5)  follows. 

APPENDIX  B 
Proof  of  Theorem  4 

The  following  nonrecursive  equation  for 
w(k)k  ‘4)  is  easily  obtained  from  Eqs.  (40),  (61), 


If  H4q  V = 0,  as  assumed,  then  N (k|k  + q)  = 0, 
which  completes  the  proof  of  the  femma. 

Proof  of  Theorem  4.  Theorem  4 is  proved  by 
repeated  application  of  our  lemma,  as  follows.  If 
HT"  = 0,  then  Nw(k|k  + 1)  = 0,  which  means  that 


1 


w(k|k  + l)  0 [sec  Eq.  (13-1) J.  U IC  O and  H*r  - 0, 
N (k|k  + 2)  = 0;  but,  since  Hr  0,  N (klk  + l)-0  as 
well;  hence  w(k|k  + 2)  = 0.  Proceeding  in  this 
manner,  we  conclude  that  if  Mr  = H)~-...  - Ht^  P 
*0.  then  Nw(k  | k + j)  * N (k|k  + j-  l)-...-N  (k  | k + 2) 

| N (k|k  + l)%V  0;  hence,  Ww(k |k  t j(  • 0.  In  s'Kort,  we 

have  shown  that,  under  the  conditions  of  Theorem 
4,  Eq.  ( 66  | is  true. 

Next,  set  q j + 1 in  Eq.  (13-9)  to  show  that 

Nw(k|k  +j  + 1)=Q  (H*V)’[HP(k  +j  + l|k  + j)H'+  R]"1 

(B  -10) 

[Since  H$^P  * 0,  N (k|k  + j+l)*0;  hence,  under  the 
conditions  of  the 'theorem, 

w(k  )k  + j + 1)  = N^(k  |k  + j + l)”S(k  + j + ljk  + j)  * £ , 

(B-U) 

which  proves  Eq.  ( 67 ),  and,  completes  the  proof. 
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Figure  1.  Example  of  reflection  coefficient 
sequence,  U(k). 
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Figure  2.  Linear  system  interpretation  for 
convolution  summation  model  of  a 
seismic  trace.  Signal  V (k)  is  given 
by  Eq.  (4). 
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